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ABSTRACT 

At ~ 400 pc, the Horsehead Nebula (B33) is the closest radiatively-sculpted pillar to the Sun, 
but the state and extent of star formation in this structure is not well understood. We present 
deep near-infrared (IRSF/SIRIUS JHKg) and mid-infrared (Spitzer /IRAC) observations of the 
Horsehead Nebula in order to characterize the star forming properties of this region and to assess 
the likelihood of triggered star formation. Infrared color-color and color-magnitude diagrams are 
used to identify young stars based on infrared excess emission and positions to the right of the 
Zero- Age Main Sequence, respectively. Of the 45 sources detected at both near- and mid-infrared 
wavelengths, three bona fide and five candidate young stars are identified in this 7'x7' region. 
Two bona fide young stars have flat infrared SEDs and are located at the western irradiated 
tip of the pillar. The spatial coincidence of the protostars at the leading edge of this elephant 
trunk is consistent with the Radiation-Driven Implosion (RDI) model of triggered star formation. 
There is no evidence, however, for sequential star formation within the immediate ~ 1.'5 (0.17 
pc) region from the cloud/H II region interface. 

Subject headings: ISM: clouds — ISM: globules — ISM: individual (Horsehead Nebula) — stars: 
formation — stars: pre-main sequence 



1. Introduction 



Star formation is a multi-faceted phenomenon that can occur spontaneously or as a reaction to external 
influences. A scenario that has gained acceptance over the last decade to account for many of the small- 
scale properties of star-forming regions involves the direct (winds and radiation) and indirect (expansion 
of H II region) triggering of star formation by massive young stars. In each case, triggered star formation 
is the result of an externally-induced compression of a molecular c loud. It has been invoked to explain 



the a pparent enhanced luminosities of pr otostars near H II regions (IMotovama et al.l 120071 : ISugitani et al 

2005), synchronized star formation 



19891 ). the age spread of young cluste rs (ILee fe Chenl l2007t ILee et al.l . 

at the peripheries of young clusters ( Smith fc BrooksF 20071). small-scale sequenti al star formation (e.g., 
Ikeda et al. 20081 : Getman et al.ll2007l Matsuvanagi et al. 200a ISugitani et al.lll995[) . and star formation in 



bright-rimmed clouds fe.g.. ISugitani fc Oguralll994l : ISugitani et al.lll99lf h 



There are tw o leading models through which triggered star formation is plausible near H II regions. 
Chen et alj (|2007f ) outline these cases and describe likely examples of this phenomenon. In the "collect and 
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collapse" model (jElmegreen fc Ladal 119771 : iHosokawa fc Inutsukal 12006) , an expanding H II region sweeps 
up material into a dense layer bordering the H II region and the parental molecular cloud. This shell of 
compressed gas a nd dust then fr agments and gravitationally collapses to form new stars. In Radiation-Driven 
Implosion (RDI: iBerto ldi 1989), the photoevaporation of a bright-rimmed cloud (BRC) by the intense UV 
radiation of nearby massive young stars creates an ionization shock-front at the surface of the cloud. The 
inward pressure initiates the formation of new cores, or the compression of pre-existing ones, that then 
collapse to form a second generation of stars. 

Evidence for triggered star formation usually relies on the spatial distribution and gradient in evolu- 
tionary phases of young stars near a BRC (RDI model) or the observation of a dense, fragmented shell of gas 
with newly formed stars surrounding an H II region (collect and collapse model) . In the case of RDI, several 
authors have discussed the difficulties of determining whether youn g stars are actual ly induced to form or 
whether they are merely being exposed by photoevaporation (e.g., iDale et al.l 120071 ). The real picture is 
probably some combination of these scenarios, where the compression of preexisting clumps may produce an 
acceleration of what would have been local spontaneous star formation. However, even before the collapse 
of a core to form a protostar, the influence of mas sive stars on dens e cores has been observed in the large 
velocity widths of cores near the M42 H II region (jlkeda et al.l 120071 ). leaving no doubt that star formation 
near massive stars is not an isolated process. 

Pillars (or elephant trunks) are commonly observed in star-forming regions and, although differing in 
shape and size, are thought to have similar formation mechanisms. As a molecular cloud is photoionized by 
a massive young star, a pre-existing dense clump will photoevaporate at a slower rate than the surrounding 
more tenuous material. A pillar will form when such a clump effectively protects the region behind it 
from photoevaporation. Although they are still being eroded away, they do so at slower rates than the 
cloud material from which they protrude. In this scenario these pillars are the astronomical analogues of 
hood oo formation s in geology. O ther theories of pillar formation include a variety of instability models 
(e-g., 



Spitzerlll954 : 



Frieman 



1954 ) and twisted magnetic flux ropes ( Carlqvist et al. 20031 ). The combination 



of stella r winds and an exp anding ionization front probably influences the formation of these structures 

These pillars point radially in the direction of the weathering source and 



as well (|Schneps et al 



1980) 



are often characterized by bright rims in Ha emission. Several studies have employed Smoothed Particle 
Hydrodynamical models in an attempt to better-understand pillar formation and evolution. In these models 
core formation and eventual core collapse occurred at the lea ding edges of pillars as a direct result of th e 
pressure caused by photoevaporation at the cloud's surface ( Gritschneder et al. 2007 : Miao et al. 2006). 
Young stars at the tips of such p illars have been ob served in over a dozen cases, including the Elephant 



Trun k Nebula (jReach et alj 120041) . the Eagle Ne bula (jLinskv et al.ll2007t ISugitani et al 



other , lesser-known pillars (jSugitani et al 
20071) . 



20021 ). the Cone Nebula ; Damn & Simon 2005 ), several pillars in 30 Doradus ( Walborn ll2002l ) , and many 



1995 



2002 



Fukuda et al 



Carlqvist et~alll2003t llJbeda fc Pellerinll2007[ I Wang et al 



At a distance of ~ 400 pc 1 Anthony- Twarog 1982h . the Horsehead Nebula (or B33: Barnard 1919 ) is 
the closest elephant trunk to the Sun and represents an ideal laboratory to study the interaction between 
massive stars and molecular clouds. It extends westward from the southern part of L1630 into the H II region 
IC 434 and points radially toward the OB system a Orionis, itself at the center of the ~ 3 Myr a Orionis 
cluster (jHernandez et al.1120071 : lOliveira et al.ll2006l ). The Horsehead has a length of roughly 0.8 pc and it is 
located about 4 pc in projection from the source of its photoerosion, a Ori A (SpT: 09.5V). This is likely 
a gross underestimation of their actual separation though; a wide range of distances to the a Ori cluster 
have been published but typically fall between 350 pc and 470 pc (jPerrvman et al.l I1997T : iHernandez et al 
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2005, de Zecuw et al. 1999), with the most recent main-sequence fitting placing this cluster at 420 ± 30 
pc (jSherrv et al.l 12008). The accurate distance to both the Horsehead and a Ori is important because the 
potenti al to trigger star form ation depends strongly upon the separation between the cloud and the massive 
star(s) (jMegeath et al ] l2002h . An historical account o f the discovery of the Horsehead Nebula and subsequent 
early plate recordings is given in lPound et al.l (|2003h . and a comprehensive account of more recent work is 
presented in Meyer et al. (in press). 

The physical characteristics of the Horsehead have been studied by m any authors. Its molecular mass 



is roughly 27 M (Pound et al 



3-7 x 10 4 cm" 3 (Phili pp et al 



2003; Philip p et alj|2006t lLada et al.lll99lik the gas density varies between 



200< 



(J in m ost of the Horsehead and rea ches a peak of 6 x 10 5 cm 3 in the 



denser of two sub-mm c ores (B33-SMM1 , Ward-T hompson et al. 2006). Its gas temperature ranges from 



10 K to roughly 20 K (jPetv et al.l 120071 : iPhilipp et al.ll2006l ). The formation and evolution of B33 and 
its characteristic substructure have been debated in the literature, but little consensus has been reached 



|Reipurth fc Bouchetllliil I Warren-Smith et al.|[l98l Neckel fc Sarcandenfl98l Iffllv-Blant et al.ll2005h 



The photodissociative properties of this region have bee n studied in depth (e.g., ICompiegne et al 



20071 : iPetv et al.ll2005t iTevssier et al.l 12004 lAbergel et alj 120031 ). but few authors have focused on its star- 
forming properties. The fir st targeted search for embedded young stars in the Horsehead was performed by 
Reipurth fc Bouchetl (j.1984). They presented shallow near- infrared images which revealed many optically- 
invisible infrared sources, one of which had a large near-infrared excess (B33-1). This same source was 
detected by IRAS (IRAS 05383-0228). Two other IRAS sources were also foun d (IRAS 05384-022 9, IRAS 
05386-0229), and it has been suggested that they might be embedded young stars ( Pound et al ]l2003h . ASCA 
observations uncovered two X-ray sources in the region ( Nakano et alJll999[ lYamauchi et al. 20001 1 but their 
status as young stars remains unclear. In this study we present near-infrared (NIR) and mid-infrared (MIR) 
imaging observations of the Horsehead Nebula in order to clarify its overall state of star formation and to 
assess the likelihood of external triggering activity. 

This paper is organized in the following manner. In § 2 we introduce the IRSF and Spitzer /IRAC 
observations, including data reduction, calibration, and analysis. The results of the NIR and MIR color- 
color diagrams and the NIR color-magnitude diagrams are described in § 3. In § 4 we discuss the current 
state of star formation within this pillar. We conclude by summarizing our results in § 5. 



2. Observations and Data Reduction 

2.1. Near-Infrared Observations 

Near-infrared observations were o btained on 13 March 2002 with the Simultane ous-3 color InfraRed 
Imager for Unbiased Survey (SIRIUS) (iNagashima et al.lll999l iNagavama et al.l 120031 ) mounted on the 1.4 
m InfraRed Survey Facility (IRSF) telescope, which is part of the South African Astronomical Observatory 
(SAAO). The SIRIUS camera allows for efficient simultaneous multi-band imaging. Focused light is split by 
two dichroic mirrors to the three J (A = 1.25 /im), H (A = 1.63 /im), and K$ (A = 2.14 /im) near-infrared 
filters. The detectors are 1024 xl024 HgCdTe (HAWAII) infrared arrays with pixel scales of $.'45 pix -1 , 
producing a 7'x7' field of view. Observations were taken in a dithered pattern to account for bad pixels and 
cosmic rays. Three sets of ten 30 s images were obtained in each filter, producing an effective integration time 
of 900 s for each coadded image. Sky frames were centered at aj20oo.o = 05:41:41, 5j2ooo.o = —02:26:30 and 
twilight flat fields were obtained at dusk. Dark current subtraction, sky su btraction, and flat field d ivision 
were performed on site using the pipeline reduction procedure outlined in INagashima et al.l (|l999f ) . The 
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reduced frames were then shifted, trimmed and median-combined using the IMALIGN and IMCOMBINE 
tasks in the NOAO/IMAGES package from the Image Reduction and Analysis Facility (IRAF) 7 . The 
combined J-, H-, and Xs-band NIR image of the Horsehead is shown in Figure [TJ 

The source detection algorithm DAOFIND in the NOAO/DIGIPHOT/APPHOT package in IRAF was 
run separately for each of the three coadded images. A visual inspection of each coadded frame was performed 
to ensure that no sources were missed and that obvious contaminants were excluded. Typical Full- Width 
at Half-Maximum (FWHM) values for the three bands were 1'.'2 and remained constant throughout the 
observations. Aperture photometry was performed using the APPHOT task in IRAF using an aperture 
radius of l'/8 with inner and outer sky annulus radii of 4'.'5 and 7'.'2, respectively. These values were chosen 
based on a series of tests using different aperture radii to optimize the signal to noise ratio of the final 
flux. A 3 a cutoff above the median value across the entire image was employed for each sky annulus to 
avoid contamination from nearby stars. Two bright sources were saturated in all three bands and one source 
corresponding to a known candidate young star was saturated in the Kg band (B33-1). 2MASS magnitudes 
for all three bands were used for these three saturated stars. 

The initial criterion for source determination was a 5 a detection in at least two bands. Later a second 
criterion that final, 2MASS-calibrated J~H and H~K$ color errors be less than 0.15 mag was employed. 
The former criterion ensured that detections over multiple bands were indeed real and was used for initial 
screening. The latter criterion was used after color transformations to the 2MASS photometric system and 
associated errors were determined. It was then applied to all sources so that the J—H and H—Ks colors 
would have reasonably small errors for placement on the NIR color-color diagram (see §2.1.1). A total of 
119 NIR sources satisfied both criteria. 

Limiting magnitudes are roughly J=19.0, H=18.0, and -Ks=17.5. As a reference, at the distance of 
B33 (roughly 400 pc), a 1 Myr brown dwarf at 0.08M Q has a Ks-b&nd magnitude of ~13 (based on the 
evolutionary tracks of Barrafe et al. 1998). This survey is therefore within reach of slightly embedded 
low-mass young stars and much more heavily embedded solar-mass young stars. 



2.1.1. Calibration Using 2MASS Photometry 



The IRSF sources were calibrated using the 2MASS catalogue ( Skrutskie et al. 200(jh . The criteria for 
selecting 2MASS point sources within the IRSF field of view included the best photometric quality flags 
(labeled "A" in the 2MASS Point Source Catalogue, with a signal to noise ratio (SNR) > 10) in at least two 
bands and a flag of at least second-best photometric quality in the third band (labeled "B," with a SNR > 
7). An additional criterion that sources have magnitudes between roughly 12 and 15 in the K$ band was 
imposed as 2MASS is most sensitive over these magnitudes. A total of 23 bright 2MASS sources were used 
to calibrate the IRSF magnitudes. 

The basic transformation equations between filter systems were used to transform colors and magni- 
tudes from the measured IRSF instrumental magnitudes to the standard 2MASS filter system (see, e.g., 
Carpenter! l200lh . The coefficients and zero-point offsets were determined by applying a linear weighted 
least-squares fit to the appropriate color and magnitude relationships. As expected, there was very little 
color dependence between the two systems. The transformation equations were applied to all instrumental 



7 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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magnitudes. SIRIUS does not record accurate sky coordinates and so all source coordinates were derived 
using the astrometry calculator Astrom 8 with the 23 brightest 2MASS stars as reference positions for the 
interpolation of coordinates. 



2.2. Mid-Infrared Observations 

The Infrared Array Camera (IRAC) on board Spitzer Space Teles cope is a four-chan nel camera that 



obtains simultaneous broad-band images in 3.6, 4.5, 5.8, and 8.0 /zm ([Fazio et al.ll2004bl ). Spitzer /IRAC 
observations of the Horsehead were obtained as part of the program titled An IRAC Survey of the LI 630 
and L1641 (Orion) Molecular Clouds (PI: S.T. Megeath; PID43). Two epochs of observations separated by 
about six months were used to construct two 12 s "high dynamic range frames." These frames consist of a 12 
s frame with 10.4 s of integration time and a 0.6 s frame with 0.4 s of integration time. The total integration 
time is 41.6 s for the long frames and 1.6 s for the short frames. The IRAC color-composite image of the 
Horsehead is shown in Figure O 

Photometry of point sources was initially obtained from the Basic Calibrated Data (BCD) from the 
Spitzer Science Center S14 data pipeline. The BCD data were mosaiced using the WCSmosaic IDL package 
developed by Robert Gutermuth. This software rejects cosmic rays by identifying deviant pixel values in 
overlapping frames, applies a frame-by-frame distortion correction, derotation, and subpixel offsetting in a 
single transformation, and matches the background by comparing signal levels in overlapping images. 

The photometry was extracted using PHOTVIS, an IDL photometry and visualization tool developed 
by Robert Gutermuth. It performs a spatial filtering of the data, computes a noise map from the filtered 
data, and then searches for point sources that exceed a source detection threshold. Aperture photometry 
was obtained for each of the sources using an aperture of 2 pixels and a sky annulus of 2-6 pixels. Each 
pixel is l'.'2x 1'.'2. The data were then bandmerged and sources were considered coincident when they were 
within l'.'O. If two sources satisfied that criterion then the closest sources were then identified as coincident. 
The zero points used to obtain the magnitudes (in DN/S per BCD image) were 19.6642, 18.9276, 16.8468, 
and 17.3909, for the 3.6, 4.5, 5.8, and 8.0 /im bands, respectively. The limiting magnitudes are [3.6] ~ 17.0, 
[4.5] ~ 16.6, [5.8] ~ 15.8, and [8.0] ~ 13.5 for errors < 0.1 mag in the first two IRAC channels and < 0.2 
mag in the last two channels. A total of 61 sources were observed with IRAC in at least one band. 



2.2.1. Rejection of Extragalactic Contaminants 



Galaxie s with high star fo rmation rates occupy the same region of the NIR color-color diagram as do 



2006h . Galaxies and young s tars also occupy a similar region of the [3.6]-[4.5] vs 
[5.8]-[8.0] color-c olor diagram ( Guterm uth et al.l 120081 ) . It is therefore imperative to remove these objects 



young stars (jGeller et al 



from the sample. iHernandez et al.l (|2007f ) describe two methods for rejecting extragalactic objects from IRAC 
data; we adopt both of these procedures in this study. Roughly half of the sour ces with [3.6] > 14.5 are 
extragalactic in nature ( Fazio et al.l 2004"al ). so we reject sources above that limit. Gutermuth et al. ( 2008h 
identify regions of the [3.6] — [5.8] vs. [4.5] — [8.0] and the [4.5] — [5.8] vs. [5.8] — [8.0] color-color diagrams that 
are dominated by polycyclic aromatic hydrocarbon-rich galaxies with active star formation. Wc eliminated 



8 Distributed through the Star link Software Collection, which is currently run by the Joint Astronomy Centre of the UK 
Science and Technology Facilities Council 
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probable contaminants by identifying sources in those regions and removing them from the sample. 



2.2.2. Merging the NIR and MIR bands 

Finally, the NIR and MIR sources were bandmerged by selecting for IRAC sources within 3'.'6 (= 
3xFWHMj_f, a , i£ i) of each IRSF point source. If more than one IRAC source satisfied that requirement then 
the nearest one was selected. The final sample used in this study consisted of 45 sources with J, H, K$ 
magnitudes and a detection in at least one IRAC band. 



3. Results 

3.1. Identification of Young Stellar Objects 

A source is considered a candidate young stellar object (YSO) if there is evidence for an excess of IR 
emission that cannot be accounted for by photospheric emission alone, indicative of a warm inner circum- 
stellar disk. Candidate YSOs are identified using the following criteria: sources must have an IR excess in 
at least two of the six NIR, MIR, and combination of NIR and MIR color-color diagrams. Alternatively, 
if the source is detected in all four IRAC channels, it is considered a candidate YSO if it has an excess in 
the MIR color-color diagram and has a spectral index > -2.56, which allows for the inclusion of transition 
or anemic disks whose spectral energy distributions only break from the photospheric emission near 8 /im. 
Furthermore, the positions of candidate YSOs in the NIR color-magnitude diagrams must be consistent with 
youth at the distance of the Horsehead (to the right of the Zero- Age Main Sequence). If a source exhibits an 
excess in a shorter-wavelength color-color diagram then it must also have an excess in the longer-wavelength 
IR color-color plots (depending on which IRAC bands it was detected in) to remain as a candidate YSO, 
as young stars with warm inner disks will have an excess at both NIR and MIR wavelengths. However, the 
converse isn't required: if the inner disk is cleared then the excess emission may only begin to appear at 
MIR wavelengths. The 3.6 to 8.0 /im spectral indices, the spectral energy distributions, and the positions 
in color-color diagrams are used together to distinguish bona fide YSOs from candidate YSOs and to assess 
infrared evolutionary classes. 



3.2. Infrared Color-Color Diagrams 

3.2.1. Near-Infrared Color- Color Diagram 



The NIR color-color diagram has long been recognized a s an important t o ol in t he identification of young 
stars with warm circumstellar disks (jLada fc Ada ms 1992). iBessell fc Brettl ([19881 ) demonstrated that main 
sequence stars and giants on the J—H vs. H—Ks diagram sit on char acteristic loci. A similar empirical 
locus of Classical T Tauri stars was later found by I Meyer et al.l (|1997T ). Interstellar reddening will affect 
the position of a star on this diagram by shifting it to redder J — H and H — Kg colors, as is indicated 
by the arrow for Ay = 5 mag in Figure [31 Sources redward of the right-most reddening envelope are good 
candidates for young stars with warm disks, although young stars with and without disks may still reside 
in the "normal" portion of this plot. We select sources in Figure [3] with H — Kg > 1 o from the right-most 
reddening envelope to be candidate young stars. There are two NIR excess sources (B33-1 and B33-28), one 
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of which was detected by 2MASS (B33-1). MIR excess sources in the [3.6]-[4.5] vs. [5.8]-[8.0] color-color 
diagram (see § 3.2.2 and Figure |U below) are overplotted as open circles. 



3.2.2. Mid-Infrared Color-Color Diagrams 

MIR color-color diagrams are useful tools to assess the IR class of young stars and can help distinguish 
betwee n highly red d ened n ormal stars and those with true IR excesses caused by emission from cirumstellar 
disks. lAllen et al.l (j2004h plotted the positions of Class I (accreting envelopes and disks) and Class II 
(accreting disks) YSO models in the [3.6] — [4.5] vs. [5.8] — [8.0] diagram over a wide range of parameters, 
including accretion rates, disk sizes, and inclinations. Unique regions were identified as being the loci of 
Class I and Class II objects. The po sitions of young stars in that diagram were subsequently shown to be in 



good agreement with these models (Mcge ath et al . 2004, 2005; Hartma nn et al.l 12005c iLuhman et al 



2006t 



Balog et al.ll2007t IWinston et al.ll2007[ ). Young stars with anemic disks or no disks at all (Weak-lined T Tauri 
stars/Class III IR sources) have similar NIR and MIR colors to main sequence dwarfs, where photospheric 
emission dominates over any slight disk emission (jHartmann et alj|2005t iHernandez et al.ll2007l ): young stars 
with highly evolved disks are therefore missed in this type of survey. 

There were 22 sources detected in all four IRAC channels (with errors < 0.1 in [3.6] and [4.5] and errors 
< 0.2 in [5.8] and [8.0]) after bandmerging with the IRSF data and after extragalactic rejection criteria were 
applied. They are plotted in the [3.6] — [4.5] vs. [5.8] — [8.0] diagram in the upper portion of Figure The 
loci for main sequence stars (or young stars w ith little excess emission) and Class II sources are plotted 
as a dotted box and a dashed box, respectively (| Allen et alj|2004t lHartmann et al J 120051 ) . Both NIR excess 
sources fall in the Class I region of the MIR color-color diagram. Their physical locations in the Horsehead 
are noteworthy: both are embedded in the western bright limb near the cloud/H II regio n interface. B33-1 
was previously identified as a young star ( Reipurth &; Bouchet 1984 ; Pound et al. 2003f ). but its IR class 
was not ascertained. Three other sources fall in the Class II region of Figure [4j B33-21, B33-25, and B33- 
31. None of the three sources exhibit a NIR excess. The [3. 6] -[4. 5] vs. [4. 5] -[5. 8] color-color diagram is 
plotted in the lower portion of FigureHJ The spread in IR classes is much smaller than that of the [3.6]-[4.5] 
vs. [5.8]-[8.0] color-color diagram, so that plot is not used to identify candidate young stars in this study. 
Nevertheless, it is instructive to note that four of the five MIR excess sources appear to have some excess 
emission in this diagram relative to the rest of the sample. The symbols are the same as the upper plot and 
all five MIR excess sources are labeled. 



3.2.3. IRSF and IRAC Combined Color- Color Diagrams 

Th e combination of NIR and MIR colors has b een shown to be useful in identifying bona fide IR excess 
sources ( Hartmann et al. 2005; Win ston et al ] l2007h . The 8.0 /mi band is particularly sensitive to PAH emis - 
sion, which has relevant emission features at 3.3, 6.2, 7.7, and 8.6 /Ltm (|Wu et al.l 120051 : iPeeters et al.l l2004). 
Color-color diagrams without the 8 /^m band are least affected by background extragalactic contaminants 
that may have been missed in previous selection criteria. The shorter wavelength IRAC channels are also 
more sensitive; consequently any young stars with warm disks that did not display a NIR excess and that 
may not have been detected in the 8 /jm channel should reveal themselves in the combined NIR and MIR 
color-color diagrams. 



We use four combined NIR/MIR color-color diagrams in this analysis: H— [3.6] vs. [3.6] — [4.5], K$ — [3.6] 
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vs. [3.6] — [4.5], J—H vs. H— [4.5], and H—K$ vs. Kg— [4.5]. The four color-color diagrams are plotted twice 
for clarity. In Figure El the sources with NIR and MIR e xcesses are emph a sized, whereas Figure El labels all 
sources with excesses in the individual NIR/MIR plots. lHartmann et al. (2005) define empirical regions of 
IR excess in the H— [3.6] vs. [3.6] — [4.5] and K$ — [3.6] vs. [3.6] — [4.5] diagrams using known young stars in 
the Taurus star-forming region. In their work there are clean breaks between Weak-lined T Tauri stars (Class 
III IR sources) and Classical T Tauri stars (Class II IR sources). These breaks are represented as dashed 
lines in the top two plots of Figures [5] and El they simultaneously display the slope of the reddening vector 
for each plot. We selected sources with [3.6] — [4.5] colors > 1 cr to the right of the dashed lines as having 
IR excesses. The bottom two plots of Figures [S] and [S] are the J—H vs. H— [4.5] and H—Ks vs. Kg— [4.5] 
color-color diagrams. The 45 sources with NIR and 4.5 /mi detections are plotted. The dashed lines mark 
the separation of sourc es with and without IR excesses and also display the slope of the reddening vector 
for each diagram ( from I Winston et al.ll2007l h Sources with H— [4.5] and Kg — [4.5] colors > 1 a to the right 
of the dashed lines (for the bottom-left and bottom-right plots, respectively, in Figures El and E]) have IR 
excesses. Twelve unique sources exhibit an excess in these diagrams, eight of which show an excess in two 
or more of the four plots. The two sources that have combined NIR and MIR excesses (B33-1 and B33-28) 
exhibit an excess in all four of the NIR/MIR color-color diagrams. 



3.3. Near-Infrared Color-Magnitude Diagrams 



In Figure [7] we present the J vs. J-K$ and K$ vs. H-K$ color-magnitude diagrams for the 45 
sources with NIR and MIR detections. The pre-main sequence evolutionary tracks of Baraffe et al. ( 19981 ) 
are overplotted. We first transformed the evolutionary tracks from the CIT system to the 2MASS system 
using the relations in ICarpenterl (|200lh and then shifted them to a distance of 400 pc. Photometry is not 
corrected for extinction. We can therefore only identify those stars which cannot be both young and at a 
distance of 400 pc by selecting against sources to the left of the Zero- Age Main Sequence (ZAMS; left-most 
isochrone at 500 Myr in Figure [7]). Two candidate YSOs lie in this region and are consequently rejected as 
being young. 



3.4. Summary of YSO Identification 

Sources are considered to be candidate YSOs if they have an IR excess in at least two color-color 
diagrams. If a source has an excess only in the MIR color-color diagram but has a spectral index indicative 
of a circumstellar disk {ajRAC > -2.56), then it too is considered a candidate YSO. Candidate YSOs with 
a shorter- wavelength color excess must also have an excess at longer wavelengths (see § 3.1 and Figures El 
Efl [6]). Bona fide YSOs are distinguished from candidate YSOs as having particularly strong IR excess 
emission and spectral energy distributions (see below) that resemble typical star + disk systems. Our final 
requirement is that candidate young stars must lie to the right of the ZAMS in the NIR color-magnitude 
diagrams (Figure [7]). A summary of the results is presented in Table CO Three of the fourteen sources with an 
excess in at least one color-color diagram are bona fide young stars, two of which (B33-1 and B33-28) have 
IR excesses and SEDs consistent with flat-spectrum objects in an evolutionary phase between Class I and 
Cl ass II. Five of the four t een ar e considered candidate young stars. We continue the nomenclature developed 



by iReipurth fc Boucheti (|1984l ) for IR sources in the Horsehcad. They began with "B33-1" and ended at 
"B33-26." Our new source identifiers begin at "B33-27" and run through "B33-33," with count increasing 
in right ascension. Only the relevant sources in this study are numbered. Photometry for the eight bona 
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fide and candidate YSOs is presented in Table 2 and their positions are overplotted on the -ff-band image of 
the Horsehead in Figure [SJ 

Spectral energy distributions (SED) for the bona fide and candidate young stars are plotted in Figure 
|U The SED of one source that did not have an IR excess but that may be the 1R counterpart to IRAS 
05384-0229 is plotted in the same figure (B33-7; see § 3.4.1 for details). The conversion between magnitudes 



Cohen et al. 


( 


2003 


) and 


Reach et al. 



(|2005l ) calibrations for 2MASS and IRAC, respectively. All SEDs are normalized to the _ff-band flux, which 
is less affected by reddening than the J-band. A blackbody spectrum with a temperature of 4000 K and 
scaled to the iJ-band flux is overplotted as a reference (dotted lines). 

The evolutionary stag e of the star plus disk /envelope syste m defines the overall shape of an SED and 
consequently the IR class ( Adams et al. 1987 ; Lada et al. 20061 ). A common MIR classification system is 
based on the value of the 3.6 to 8.0 /im spectral index, defined as a — d\og XF\/d\og A. In this study we 
generally use the following classification scheme, although we also take into account the positions in color- 
color diagrams to assign an IR class (Tabled]): IRAC Class I sources have olirac > 0-3, IRAC flat-spectrum 
sources have —0.3 < olirac < 0-3, IRAC Class II sources have —1.8 < olirac < —0.3, sources wit h transition 
anemic") disks have —2.56 < oirac < —1-8, and IRAC Class III sources have olirac < —2.56 (jLada et al 



2006). Spectral indices are calculated by fitting a line to the four IRAC fluxes in log(Ai<A) vs. log(A) space. 
The slope of the line is the spectral index. Values of the spectral indices for c andidate young st ars are listed 
in Tabic [T] Note that an M-dwarf photosphere has a spectral index of -2.66 ( Lada et al. 20061 ). 



3.4-1- Counterparts to IRAS and ASCA Sources 

Protostars have warm disk/envelope systems that reprocess stellar light into longer wavelength radiation. 
They are therefore often identified in FIR surveys of star- forming regions. More evolved young stars both 
with and without disks exhibit X-ray emission likely caused by intense magnetic activity generated by a 
strong dynamo effect. Large flaring events are thought to be the origin of the typical X-ray variability of 
Classical T Tauri stars. The detection of both FIR IRAS sources and ASCA X-ray sources in the vicinity 
of the Horsehead Nebula supports the notion that there is ongoing star formation in this region. It has 
been suspected that the IRAS and ASCA sources are young stars, but large positional uncertainties in those 
data have prevented the identification of most of their optical/IR counterparts. Three IRAS sources were 
detected in the Horsehead : IRAS 05383-0228, IRAS 05384-0229, and IRAS 05386-0229 (photometry is 
given m bound et al1l2003h . B33-1 has been identified as the optical/IR counterpart of IRAS 05383-0228. 



IRAS 05384-0229 was only detected in the 100 /im band and likely originates from infrared cirrus: all three 
sources have contamination flags that indicate a high probability of being caused by cirrus emission. These 
contamination flags are based on a large number of nearby 100 //m-only sources, the presence of large-scale 
structure in the 100 /jm emission map, and a high surface brightness at 100 ^m in the region. IRAS 05386- 
0229 was only detected in the 12 and 25 /im bands of IRAS. It too has high cirrus contamination flags. 
Nevertheless, we will compare the positions of the three IRAS sources with the locations of the candidate 
young stars identified in this work to search for possible YSO counterparts. Two X-ra y sources in thi s regio n 



were detected by AS CA by two different gro ups: [NY99] C-18 and [NY99] C-20 by iNakano et all (|l999f ). 



and [YKK2000] A4 bv lYamauchi et all (|2000f ). [NY99] C-18 was detected by the Gas Imaging Spectrometer 
(GIS) and the Solid-State Imaging Spectrometer onboard ASCA. [NY99] C-20 was detected by the CIS 
only, as was [YKK2000] A4. It is likely that [NY99] C-18 and [YKK2000] A4 are the same source as their 



-10- 



separation is only 18'.'5 while their positional uncertainties are both about 1'. Yamauchi et al. ( 2000h suggest 
B33-3, B33-4, B33-5, B33-9, and B33-10 as possible counterparts to [YKK2000] A4. 

In Figure [TU] we plot the positional uncertainties of the IRAS and ASCA observations. The upper 
plot shows the positions of the bona fide and candidate YSOs, while the lower plot displays all the NIR 
and MIR source detections. In the lower plot, IRSF and IRAC detections before merging the two data sets 
and before removing extragalactic contaminants are displayed. We required simply that the IRAC sources 
have errors < 0.1 in channels 1 and 2 and errors < 0.2 in channels 3 and 4. The IRSF sources with color 



errors < 0.15 are plotted. We used a 1' error circle for the ASCA observations (jYamauchi et al.l l2000). We 
searched the archival ROSAT database and found that two pointed observations were made centered near the 
Horsehead with the Position Sensitive Proportional Counter (PSPC) and the High Resolution Imager (HRI) 
detectors. The observation sequence for the PSPC frame is RP900189 and for the HRI frame is RH201148 
(PI: J.H.M.M. Schmitt); observations were made between 19 and 20 Sept 1991 and between 12 and 14 Sept 
1992 for the PSPC and HRI frames, respectively. Two X-ray sources are detected nearly on-axis in the 
PSPC observation (RX J054058-0225.5 and RX J054106.8-022349) and only one was detected with the HRI 
(RX J054106. 8-022349). We use the IAU naming convention for both sources. This non-detection of RX 
J054058-0225.5 in the HRI observation suggests that it is variable in X-rays. We determined the centroid of 
each source and plotted their positional unce rtainties in Figure [10"| The centroids are g iven in Table [3l We 
used e rror circles of 30" for the PSPC frame (jBarbera et al.lbooj iBocchino et al.ll200ll ) and 6" for the HRI 
frame ( Bocchino et al.l l2001). 



We invoke Ockham's razor and assume that the [YKK2000] A4, the NY[99] C-18, and RXJ054058- 
0225.5 are the same source, as are [NY99] C-20 and RXJ054106.8-022349. We identify B33-10 as the likely 
counterpart of the RXJ054058-02225.5 group of X-ray sources (see Figure [TU1 bottom), which sits neatly in 
all three error circles. Similarly, we identify B33-32 as the likely counterpart of the other group of X-ray 
sources centered around RXJ054106. 8-022349. It is the nearest IR source to the centroid of RXJ054106.8- 
022349 out of the two IR sources detected by both IRSF and IRAC within the RXJ054106.8-022349 error 
circle, the other source being B33-9 (not labeled). B33-1 is re-identified as the likely counterpart to IRAS 
05383-0228. No candidate YSOs are located in the error ellipses of the other two IRAS sources, but B33-7 
lies near the centroid of the IRAS 05384-0229 detection. While this IRAS source is likely caused by cirrus 
emission and is probably not a real point source, we simply note that B33-7 lies near the center of the IRAS 
error ellipse. No counterparts arc identified for IRAS 05386-0229. These results are tabulated in Table [3] 



4. Discussion 

We identify two bona fide YSOs (B33-1 and B33-28) at the western limb of the Horsehead Nebula 
as flat-spectrum protostars based on their positions in color-color diagrams, their spectral indices, and the 
overall shape of their SEDs. That their SEDs are so similar argues that they are in the same evolutionary 
phase, and consequently that they formed at nearly the same time. Their apparent separation is 26", or ~ 
1040 AU (0.05 pc) at 400 pc. B33-1 appears to be emerging from an optically-visible cavity (see Figure 
5 m bound et alll2003l ). Immediately west of the cavity is a small filamentary structure that is barely 



detected in our NIR images. It is possible that this thin filament is an an outflow from B33-1. B33-28 is 
optically-invisible and so is either embedded or is emerging on the other side of this nearly edge-on pillar. 

Ward-Thompson et al. I (|2006h recently presented observations of the Horsehead in 450 /jm and 850 /jm 



in which they discussed the stability of two dense condensations found in the neck and the head regions 
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(B33-SMM1 and B33-SMM2). They calculate the virial equilibrium of both clumps and conclude that B33- 
SMM1 has been affected by the ionizing radiation of a Ori and may undergo triggered star formation, while 
B33-SMM2 is likely a pre-existing clump in gravitational equilibrium and may be a protostellar core. In 
Figure [TT] we overplot the locations of both clumps along with the bona fide and candidate young stars from 
this study to search for any spatial coincidence. The bes t-fitting ellipses that encomp ass the 3 a sub-mm 
contours are drawn as solid lines (parameters are given in I Ward-Thompson et al.ll2006l ). B33-1 and B33-28 
are located at the periphery of the B33-SMM1 condensation. Likewise, B33-31 is not centered on the B33- 
SMM2 clump, although it does border that region. B33-31 is apparently an extremely embedded source 
given its highly reddened positions in all of the color-color diagrams. We estimate a reddening of Ay ~ 22 
mag from its location in the NIR color-color diagram down to the locus of Classical T Tauri stars. That the 
centroids of B33-31 and B33-SMM2 do not coincide bolsters the notion that B33-SMM2 is a pre-stellar core 
that has yet to undergo star formation. 

The overall picture of star formation in the vicinity of the Horsehead is not straightforward. It appears 
that several generations of young stars occupy the line-of-sight path to this pillar, but it remains unclear 
which sources are actually associated with the Horsehead (besides the protostars). One way to better 
understand the spatial proximity of candidate YSOs to the Horsehead is to look at their reddening in the 
NIR color-color diagram. If one assumes that B33-21 and B33-25 should lie near the Classical T Tauri star 
locus in Figure [31 then both sources appear to be in the foreground of the cloud as there is virtually no 
indication of reddening. The X-ray source B33-10 is slightly reddened (Ay ~ 4 mag) and may be partially 
embedded or may be a background source seen through tenuous material at the "nose" of the Horsehead. 
In the latter scenario it is possible that B33-10 is a p art of the extended a O ri cluster, as members have 
been found almost all the way to the Horsehead itself ([Hernandez et al.l 120071 ). In fact, it could be argued 
that the newest members of the a Ori cluster are just emerging at the tip of the Horsehead. This becomes 
somewhat arbitrary though as all these disparate populations eventually blend into each other. The other 
X-ray source B33-32 sits directly on the mid-M portion of the main-sequence locus in Figure [3] This raises 
the possibility that B33-32 might be a foreground M-dwarf, as magnetic M-dwarfs are strong X-ray emitters. 
Follow-up spectroscopic observations are required to fully understand most of the young stellar population 
in this region. 

Unfortunately, this study sheds little light on the origin of the unique structure of the Horsehead. 
Nevertheless, we are able to state that there is no evidence of an IR excess in the previo usly suspected young 



star B33-6, which it was suggested had formed the "jaw" cavity through an outflow (jReipurth k. Bouchet 
1984J). Neither do we find evidence for small-scale sequential star formation within ~ 1.'5 (0.17 pc) of the tip 



of the Horsehead, although we do find simultaneous star formation at its leading edge. It is probably not 
possible to conclusively state that the Horsehead is a site of triggered star formation as there are very few non- 
degenerate predictions of spontaneous and triggered star formation theories. Nevertheless, the observation 
that most relatively massive irradiated elephant trunks only show young stars at their tips suggests that 
what we are observing is indeed the result of a triggering proc ess. Noting that the sum-mm clump B33- 
SMM1 was probably formed from the direct influence of a Ori ( Ward-Thompson et al. 2006) and that the 
two simultaneously-formed protostars positionally coincide with the leading edge of this pillar suggests that 
these protostars were p robably formed as a result of the influence of a Ori. The formation timescale of the 
Horsehead is ~ 10° y rs (bound et al.ll2 0031 which is similar to the estimated ages of flat-spectrum protostars 
( Kenvon et al. 1990 : Greene fc Mever 1995h . These similar timescales are consistent with the picture of 
B33-1 and B33-28 forming when the pillar itself began to take shape. Our observations closely resemble 
the simulations of iGritschneder et alJ (j2007h . who show that ionizing radiation can lead to the formation of 
pillars similar in size to the Horsehead, and that core formation and collapse at the tips of such pillars can 
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occur. 



5. Conclusion 

We present near-infrared (IRSF/SIRIUS) and mid-infrared (Spitzer/IRAC) observations of the Horse- 
head Nebula in order to study the state and extent of star formation in this nearby pillar. We use six near- to 
mid-infrared color-color diagrams to select candidate young stars based on their infrared excesses. We also 
ensure that their positions in near-infrared color-magnitude diagrams are consistent with youth. In total we 
find two flat-spectrum protostars (B33-1, B33-28), a Class II YSO, and five candidate young stars. Spectral 
energy distributions indicate thick disk/envelope systems for the protostars and only modest excesses in 
most of the candidate young stars. We identify the infrared counterparts of two X-ray sources in the region 
and the counterpart to one of three IRAS sources, the other two likely being caused by IR cirrus emission. 

Star formation in the Horsehead is ongoing and is likely triggered by the influence of the nearby a 
Orionis OB system. We find no evidence for sequential star formation in this region, but the two proto- 
stars are in the same evolutionary phase and were probably formed from the collapse of a core that has 
since been photoevaporated away. The formation timescale of the Horsehead is also similar to the esti- 
mated ages of protostars. These observations are consistent with the Radiation-Driven Implosion model of 
triggered star formation, although it is not possible to tell whether star formation in this region is simply 
an accelerated version of what would have occurred naturally, or whether star formation is taking place 
where it otherwise not have spontaneously done so. However, one available clue is the observation of an 
externally-compr essed sub-mm clump at the w estern limb of the Horsehead (B33-SMM1), whose properties 



are discussed bv IWard-Thompson et al.l (|2006l ). This clump is out of virial equilibrium and, as such, may 
have been compressed by the influence of a Ori. That such an influence is evident immediately behind the 
locations of these protostars may suggest that they were formed in a similar fashion. 

The candidate young stars surrounding the Horsehead range from foreground to partly embedded to 
fully embedded objects. Their association with the Horsehead remains unclear. Follow-up spectroscopy and 
longer-wavelength photometry will be useful to confirm their status as young stars. 
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Fig. 1. — Three-color composite near-infrared image of the Horsehead Nebula. Blue, green, and red represent 
the J-, H-, and i^g-band IRSF/SIRIUS images. The images have been Gaussian smoothed for better 
rendering. 
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Fig. 2. — Spitzer IRAC image of the Horsehead Nebula. The image is a three color composite constructed 
from the 3.6 /im (blue), 4.5 /im (green) and 8.0 /im (red) images. The extensive nebular emission is due to 
UV heated hydrocarbons. 
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Fig. 3. — J — H vs. H — K$ NIR color-color diagram. Sources > 1 a to the right of the right-most reddening 
envelope [dotted lines) have an excess of NIR emission, likely caused by a warm circumstellar disk. They 
are plotted as filled black circles. Solid lines represent the empirical locus of main sequence stars and giants 
as determined by Bessell & Brett (1988). An analagous T Tauri Star locus was found by Meyer, Calvet, & 
Hillenbrand (1997) and is plotted as a dashed line. MIR excess sources from the 4-channel IRAC color-color 
digram (Figured! top) are overplotted as open circles and are labe led. All sources w ithout a NIR excess are 
plotted as filled gray circles. The reddening vector is derived from Tokunaga ( 2000l ). 
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Fig. 4. — [3.6]-[4.5] vs. [5.8]-[8.0] color-color diagram for 22 sources detected in all four IRAC channels 
(top). The five sources that exhibit MIR excesses are labeled (open circles). Three fall in the Class II region 
and two fall in the Class I region (from Allen et al. 2004, Hartmann et al. 2005, and Luhman et al. 2006). 
The two NIR excess sources from Figure [3] are plotted as filled black circles. Those with no MIR excess 
are plotted as filled gray triangles. Sources detected in the first three channels of IRAC arc plotted in the 
[3.6]-[4.5] vs. [4.5]-[8.0] diagram (bottom). The spread between IR evolutionary classes is smaller than that 
for the 4-channel color-color diagram, so this plot is not used to select YSO candidates; it is nevertheless 
instructive to show the similar trend of a red [4.5]-[5.8] color for candidate YSOs. Symbols are the same as 
in the upper diagram. 
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Fig. 5. — NIR/MIR color color diagrams showing the positions of NIR excess sources from Figure |3] (filled 
black circles) and MIR sources from Figure [4] (open circles). All sources not exhibiting an IR excess are 
plotted as filled gray circles. The dashed lines represent the break between normal and IR excess stars 
(Hartmann et al. 2005; Winston et al. 2007). The slope of the dashed line is the direction of the reddening 
vector for each diagram. 
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Fig. 6. — Similar to Figure [Bj but emphasizing the IR excess sources in each diagram. Sources exhibiting an 
IR excess are plotted as open triangles and are labeled; all others are represented as filled squares. There 
are 12 unique sources that exhibit an excess in at least one diagram. 
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Fig. 7. — NIR color-magnitude diagrams showing the positions of the 45 sources with NIR and MIR detec- 
tions. The pre-main sequence evolutionary tracks are from Baraffe et al. (1998). Isochrones are plotted as 
solid lines for 1, 4, 10, 40, and 500 Myr (ZAMS) from upper right to lower left. Iso-mass tracks are plotted 
as dotted lines for 0.08, 0.2, 0.4, 0.8, and 1.4 M Q (from lower right to upper left). The tracks were converted 
from the CIT photometric system to the 2MASS system using the relations in Carpenter (2001). They were 
then shifted to a distance of 400 pc. Sources with an IR excess in at least one color-color diagram are plotted 
as open squares; those with no excess are plotted as filled squares. The two sources to the left of the ZAMS 
are inconsistent with being young and at a distance of the Horsehead. They are therefore rejected as being 
candidate YSOs. Many of the IR excess sources show large amounts of reddening. The reddening vector is 
derived from Tokunaga (2000). 
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Fig. 8. — Locations of bona fide YSOs (circles) and candidate YSOs (diamonds) overplotted on the ff-band 
image of the Horsehead. The two fiat-spectrum sources are located at the western interface of the cloud/H 
II region (B33-1 and B33-28). 
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Fig. 9.— SEDs of bona fide YSOs ("YSO"), candidate YSOs ("CYSO"), X-ray sources ("XRAY"), and/or 
far-infrared sources ("FIR"). IRSF, IRAC, and IRAS data are labelled as filled triangles, filled circles, and 
filled diamonds, respectively. Magnitudes are con verted to fluxes usin g the absolute ca libration of Rieke et al 



(2008), who modified previous calibrations from lCohen et al.l | |2003h for 2MASS and lReach et alj |2005l ) for 
IRAC. The zero-point offsets used were 1617, 1073, 675, 280.9, 179.7, 115.0, and 65.09 Jy for the 1.25, 1.65, 
2.15, 3.6, 4.5, 5.8, and 8.0 /im bands. SEDs were normalized to the ff-band. Dotted lines display a 4000 K 
blackbody, also normalized to the f/-band data point. Note the odd 100 /im point of B33-7 (which does not 
show an excess in any color-color diagram); this point is probably caused by coincident IR cirrus emission. 
The IRAS data for this source was flagged as having a strong likelihood of contamination. 



- 26 - 



-2.40 



-2.42 - 



-2.44 



-2.46 



-2.48 - 



J5-2.50 



-i 1 1 1 1 1— I 1 1 1 1 r^-i 1 1 1 1 1 1 1 1 1 1- 



[NY99] C-20 



/ 

y' 

RX J054106.8-022349_ 

O B33-25 



[YKK2000] A4. 



[NY99] C-18 ( | -1 ♦B33 X -io I / 




O B33-22 

IRAS 05386 

O B33-33 <g, B33-2I 



RX J054058-0225.5 

♦ B33-31 , RAS 05384 




-2.40 " 



-2.42 



-2.44 



-2.46 



-2.48 



-2.50 k 



-i 1 1 1 1 1 1 1 1 1 1 r^-i 1 1 1 1 I A | 1 1 1 T 

H < ' A ^ ,A A 3 A 



' 
\v /A 



/ 



I ' , A H s \ 



-.A 



I 



) A A fj A 
B33-10--\ V . ■ , I / □ 

'• y A 



□ ^ ^ A 

333-1 ^ 



A AA 



85.30 85.28 85.26 85.24 85.22 85.20 
«2000.0 ( de 9) 



Fig. 10. — Identification of X-ray and far-IR counterparts. In the upper panel we plot the error circles 
(ellipses for IRAS sources) for the two X-ray sources detected by both ASCA and ROSAT (from archival 
observations) together with the three IRAS sources in the field. Overplotted as filled diamonds are the bona 
fide and candidate young stars identified in this survey. In the lower panel all IRSF sources with color errors 
< 0.15 mag (open triangles) and IRAC sources with errors < 0.1 mag in the first two bands and errors < 0.2 
mag in the last two bands (open squares) are overplotted in gray. We identify possible IR counterparts to 
both X-ray sources and to two IRAS sources. However, one IRAS source for which we identify a counterpart 
(IRAS 05384-0229; B33-7) was only detected in the 100 ^m band and has a high probability flag for being 
cirrus emission. The abnormal SED of B33-7 with the IRAS 100 /im flux supports the idea that this is likely 
caused by IR cirrus. 
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Fig. 11. — A comparison between the two sub-mm clumps identified by Ward-Thompson et al. (2006) and 
the bona fide and candidate young stars identified in this survey (filled diamonds). The solid lines show 
the best-fitting ellipses to the 3 a sub-mm contours of the B33-SMM1 and B33-SMM2 clumps. B33-1 and 
B33-28 are near the periphery of the centroid of B33-SMM1, as is B33-31 for B33-SMM2. We expect sources 
that may be the products of star formation in these clumps to be in near the clumps' centroids; we do not 
identify any young stars that are associated with the clumps. 
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Tabic 1. Summary of Color-Color and Color-Magnitude Diagrams 



CKJ2000.0 (5,72000-0 YSO IR 



ID 


(h m s) 


(o ■ ») 


CC1 


CC2 


CC3 


CC4 


CC5 


CC6 


CM1 


CM2 


aiRAC 


Status 


Class 


B33-27 


05 40 48.20 


-02 25 37.8 


N 


N 


Y 


Y 


N 


N 


N 


N 


-2.65 


N 




B33-1" 


05 40 51.70 


-02 26 48.7 


Y 


Y 


Y 


Y 


Y 


Y 


Y 


Y 


0.20 


YSO 


FS 


B33-28 


05 40 52.40 


-02 27 12.7 


Y 


Y 


Y 


Y 


Y 


Y 


Y 


Y 


0.21 


YSO 


FS 


B33-29 


05 40 57.60 


-02 24 09.3 


N 




N 


N 


Y 


Y 


Y 


Y 




N 




B33-10 a 


05 40 58.68 


-02 25 26.8 


N 


N 


N 


N 


N 


Y 


Y 


Y 


-2.48 


CYSO 


III 


B33-13" 


05 41 01.14 


-02 28 12.0 


N 




N 


Y 


N 


N 


Y 


Y 




N 




B33-20 a 


05 41 09.97 


-02 29 11.6 


N 




N 


N 


N 


Y 


Y 


N 




N 




B33-30 


05 41 02.09 


-02 23 44.1 


N 




Y 


Y 


N 


N 


Y 


Y 




N 




B33-31 


05 41 05.53 


-02 27 17.9 


N 


Y 


N 


N 


N 


N 


Y 


Y 


-0.84 


YSO 


II 


B33-21 a 


05 41 10.31 


-02 29 03.9 


N 


Y 


Y 


Y 


N 


N 


Y 


Y 


-1.97 


CYSO 


II 


B33-22 a 


05 41 12.23 


-02 27 34.3 


N 


N 


Y 


Y 


N 


N 


Y 


Y 


-2.31 


CYSO 


TD 


B33-23 a 


05 41 13.04 


-02 29 12.2 


N 


N 


N 


Y 


N 


N 


Y 


Y 


-2.80 


N 




B33-25 a 


05 41 14.06 


-02 25 05.1 


N 


Y 


N 


N 


N 


N 


Y 


Y 


-1.62 


CYSO 


TD 


B33-33 


05 41 14.12 


-02 29 01.6 


N 








Y 


Y 


Y 


Y 




CYSO 





Note. — All sources with an IR excess in at least one color-color diagram are listed. The color-color diagram numbering system 
is as follows: CC1 = J-H vs. H-K s , CC2=[3.6]-[4.5] vs. [5.8]-[8.0], CC3=//-[3.6] vs. [3.6]-[4.5], CC4=K s -[3.6] vs. [3.6]-[4.5], 
CC5— J — H vs. H— [4.5], and CC6— H — Ks vs. Kg— [4.5]. A "Y" for the color-color diagrams indicates that the source exhibits an 
excess in that plot. Color-magnitude diagrams are as follows: CM1-J vs. J—Ks, CM2—Ks vs. H-Ks- A "Y" in the color-magnitude 
columns means that the position of that source is consistent with being young and at a distance of the Horschcad (i.e. to the right of 
the ZAMS). Spectral indices in Column 12 are simply the slope of the best-fitting line to the 4-channel IRAC data in the log(AF^) 
vs. log(A) space. Column 13 indicates the status of the sources as YSOs. Bona fide YSOs are labelled as "YSO" while candidate 
YSOs have a "CYSO" label. Sources that are not YSOs arc labeled "N" . The infrared class is listed in Column 14 for flat-spectrum 
sources ("FS"), Class II sources ("II"), transition disk objects ("TD"), and Class III sources ("III"), which were determined using a 
combination of SEDs, positions in color-color diagrams, and o^rac indices.) 

a Sourcc numbers from Rcipurth &c Bally (1984) 
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Tabic 2. Infrared Photometry 



ID J (crj) J-H(ctj-h) K s (<tk s ) H - K s (<jh-k s ) [3.6] (<t 3 .6) [4.5] (0-4.5) [5.8] (cts.s) [8.0] (crg.o) 

Bona Fide YSOs 

B33-l a 12.25 (0.04) 1.52 (0.05) 9.71 (0.03) 1.03 (0.04) 8.362 (0.002) 7.586 (0.002) 6.824 (0.002) 5.717 (0.002) 

B33-28 15.79 (0.04) 1.10 (0.05) 13.53 (0.06) 1.08 (0.03) 11.407 (0.018) 10.781 (0.007) 9.947 (0.056) 8.792 (0.098) 

B33-31 18.98 (0.14) 2.97 (0.15) 14.36 (0.08) 1.57 (0.05) 13.180 (0.010) 12.908 (0.009) 12.475 (0.076) 11.464 (0.159) 

Candidate YSOs 



B33-10 


12.30 (0.04) 


1.42 (0.05) 


10.52 (0.02) 


0.38 (0.03) 


9.583 (0.002) 


9.484 (0.002) 


9.358 (0.005) 


9.294 (0.029) 


B33-21 


13.85 (0.02) 


0.58 (0.04) 


12.98 (0.02) 


0.29 (0.03) 


12.659 (0.005) 


12.499 (0.005) 


12.385 (0.031) 


11.898 (0.073) 


B33-22 


13.66 (0.04) 


1.59 (0.06) 


11.37 (0.04) 


0.68 (0.02) 


11.096 (0.005) 


10.852 (0.004) 


10.799 (0.036) 


10.616 (0.153) 


B33-25 


12.80 (0.02) 


0.69 (0.04) 


11.92 (0.02) 


0.20 (0.04) 


11.791 (0.007) 


11.796 (0.006) 


11.425 (0.064) 


10.798 (0.152) 


B33-33 


17.03 (0.05) 


0.97 (0.06) 


15.71 (0.02) 


0.35 (0.06) 




14.528 (0.041) 






X-ray Counterparts 11 


B33-32 


13.04 (0.02) 


0.61 (0.04) 


12.17 (0.02) 


0.26 (0.03) 


11.925 (0.007) 


11.918 (0.006) 


11.808 (0.068) 





a Saturatcd in the NIR images. JHKs photometry is from 2MASS. 

b B33-10 
this table. 
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Table 3. Likely Counterparts to Far-IR/X-Ray Sources in the Horsehead 



IR 

Counterpart 


a J2000.0 

(h m s) 


<5j2000.0 
(0 ' ») 


aiRAC 


YSO 
Status 


Source ID 


Type of 
Source 


B33-l a 


05 40 48.20 


-02 25 37.8 


0.20 


YSO 


IRAS 05383-0228 


FIR 


B33 _7 ? a,6 


05 40 57.26 


-02 27 44.0 


-2.84 


N 


IRAS 05384-0229 


FIR 


B33-10 a 


05 40 58.68 


-02 25 26.7 


-2.48 


CYSO 


RX J054058-0225.5 


X 












[NY99] C-18 


X 












[YKK2000] A4 


X 


B33-32 


05 41 06.89 


-02 23 49.6 




X 


RXJ054106. 8-022349 


X 












[NY99] C-20 


X 



Note. — Probable infrared counterparts to far-infrared (FIR) and X-ray sources. Labels for YSO 
status arc the same as for Table [T] 

a Sourcc numbers from Rcipurth & Bally (1984) 

b B33-7 lies near the centroid of IRAS 05384—0229, but that IRAS source is probably caused by cirrus 
emission. See § 3.4.1 for details 



